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Abstract 

A model for diffusion on a cubic lattice with a random distribution of 
traps is developed. The traps are redistributed at certain time intervals. 
Such models are useful for describing systems showing dynamic disor- 
der, such as ion-conducting polymers . In the present model the traps 
are infinite, unlike an earlier version with finite traps, this model has a 
percolation threshold. For the infinite trap version a simple analytical 
calculation is possible and the results agree qualitatively with simulation. 

PACS nos : 05.40.+j, 66.30.Dn, 61.43.-j 

1 Introduction 

Diffusion in a disordered system is a well-studied problem |], . Different 
approaches have been taken - such as models which consider inter-site 
barriers of varying heights, or models which picture sites for diffusing 
particles as wells with varying depths, also combinations of barrier and 
well models 0. Problems where the diffusing particle may encounter 
randomly distributed traps are also studied ||, [5], ||. However, one class 
of problems involve dynamic disorder, i.e. wells, barriers or traps whose 
distribution in space changes with time. In other words the location 
of the traps get redistributed at certain time intervals. Such models 
describe systems such as ion-conducting polymers, or glasses above the 
glass transition temperature |J. 

Models with dynamic disorder have also been studied for some time. 
The 'dynamic bond percolation model' (DBPM)^, [| by the Northwest- 
ern University group, and its many variations study one aspect of the 
problem quite extensively, effective medium approaches have also been 
used [TTJ, [TTJ. 



A different formulation of the dynamic disorder problem, considering 
'well' model, rather than a 'barrier' model was suggested by Bhat- 



tacharyya and Tarafdar |12|]. Mandal et al [Q applied model [|12j to ex- 



plain experimental results on the PEO-NH4CIO4 system [13]. The model 
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1^1 considers a random distribution of two types of sites on a square 



lattice, both having finite dwell times, i.e. non-zero jump rates. It is not 
a 'percolation' model in the sense that it never has zero diffusivity. The 
diffusion coefficient for different ratios of the two types of sites, and for 
different rearrangement times was calculated by computer simulation. 

In the present work we take the rearranging lattice model with finite 
traps [ 12 1 to the limit where one type of site is an infinite trap. In this 
limit the diffusion coefficient as function of the trap concentration c and 
rearrangement time r r are calculated in three-dimensions. The results are 
compared with the finite trap model as well as DBPM, and are supported 
by computer simulations. 



2 The rearranging trap model 

We consider a cubic lattice with a fraction c of sites occupied by traps. A 
particular configuration of traps remains constant from t = to a time 
interval t = r r . After that, the traps get rearranged, though the trap 
concentration c is constant. If Nq walkers started a random walk on the 
lattice from different sites at t — 0, during time r r a certain number get 
trapped. The total square distance that would have been travelled by 
the No walkers in the absence of traps is 

(r 2 ) tot = 2.d.N .D r r (1) 

If there are traps this value gets reduced to say (r tr 2 )tot- Here D Q is the 
diffusion coefficient with no traps, and d is the dimension of the system. 

After the interval r r , since traps are rearranged, previously trapped 
walkers may be set free and resume their walk. In the next interval from 
t = T r to t — 2r r they again cover a total distance(r 4r 2 ) to i. Fig 1 shows 
how the square distance covered increases with time in our simulation de- 
scribed in section 3. So for a time interval >> r r , the diffusion coefficient 
on a rearranging lattice is given by 

D tr = (r tr 2 ) tot /2.d.N .r r (2) 

This argument is the same as in M. The problem is thus reduced to 
calculating (r tr 2 ) io t.With traps having infinite depth (i.e. infinite dwell 
time), this is estimated as follows, without involving any complicated 
mathematics. 

2.1 Calculating the diffusion coefficient D(c,r r ) 

Suppose Nq particles start a random walk at t = on a cubic lattice 
occupied randomly by c fraction of traps. The remaining (1 — c) sites are 
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'allowed' sites. If a particle encounters a trap it gets stuck until the next 
rearrangement time , (r tr 2 ) is calculated as the mean distance travelled 
by each of the N particles (r 2 tot )/N . If there were no traps (r 2 ) tot 
would be given by equation 1. This is the total distance travelled by all 
the particles in time r r . In presence of traps more and more particles get 
stuck and the total square distance travelled does not reach the value in 
eqn (1), but falls short by an amount Ar. 

Ar = 2.d.N .D .r r - {r 2 tr ) tot (3) 

knowing Ar one can calculate D(c, r r ) from the relation 

D(c,r r ) = (r tr ) 2 tot /2.d.N Q .T r (4) 



2.2 Calculation of Ar : the trapping law 

To calculate Ar we assume a simple law for trapping. 
We assume that the number of particles dN trapped in the time interval 
t to t + dt is proportional to N - the number of particles at time t, c - 
the trap concentration and dt. So, 

dN = -a.c.N.dt (5) 

a ~ dS(t)/dt, where S(t) is the number of distinct sites visited in time t 
E| . For three or higher dimensions 

S{t) oc t 

[y| , so a is a constant depending on the dimension of the system but 



independent of t. We have thus the simple trapping rule 

N = N .exp(-a.c.t) (6) 

It may be mentioned here, that problems of absorption by traps, is 
well studied, and exact analytical results have been worked out for the 
asymptotic limit t — > oo JTB], [T7|. In the asymptotic limit the number of 
surviving particles at time t is given by 

N = N .exp(-a.c 2 / {d+2 \t d ^) (7) 

where the exponents of c and t as well as a are dimension dependent. 
However this result is valid only when the probability of survival becomes 
very small ~ 1(T 13 HSU. In the regime we are interested in the simple 
expression eq(6) agrees much better with simulation results than eq (7) as 
we show in fig 2 for trapping in three-dimensions. We may now proceed 
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to calculate D. The particles trapped at time t to t = t + dt contribute 
an amount 5r less to (rf r ) tot than they would if not trapped. 

Sr = 2.d.(r r - t) (8) 

Henceforth we consider the case d = 3 only. 

Integrating this over the interval r r for all particles trapped during 
this time, we have 

Ar = / 5rdN = 6 / a.c.N .exp(—a.c.t)dt.{r r — t).Do (9) 
Jt=o Jo 

or, 

= 6.D .N (r r + exp(—a.c.r r )/a.c — 1/o.c) (10) 

Since 

(rl) tot = Q.N .D r r - Ar (11) 

the diffusion coefficient in the presence of traps with a rearrangement 
time r r is 

D = (D /acr r )(l - e' ac ^)(l - c) (12) 

the factor (1 — c) accounts for the particles which are trapped at the 
outset i.e. at t — 0. 



2.3 The limiting cases 

For T r — > that is for very rapid rearrangement 

D/D = 1 - c (13) 

as expected. For r r — > oo that is when the system is effectively quenched, 
with no rearrangement 

D/D = 1 (14) 

for c = and 

D/D = (15) 

for c > 0. So there is a percolation threshold at c = 0. Any nonzero 
trap concentration, however small gives D = in the limit r r — > 0. Fig 
3 shows D vs c for different r r . The value of a is taken as a = .68, as 
given by our simulation results for the survival of random walkers in a 
three-dimensional lattice with traps. 
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3 Simulating diffusion on the rearranging 
trap model 

We have simulated the random walk on a three-dimensional rearranging 
lattice, and calculated the diffusion coefficient for different c and r r . The 
algorithm is somewhat similar to the finite trap model in fP2|. 



We work on an effectively infinite lattice, thus avoiding finite size ef- 
fects. A random walker starts to walk on a cubic lattice of unit spacing 
with a concentration c of traps. The location of the traps are not pre- 
assigned, the walker decides whether the current site is a trap from a 
random number, as it goes on. However, once assigned, the trap location 
remains fixed for time r r until the next rearrangement. Each trap is of 
infinite depth, so the walker encountering a trap has to remain there until 
the next rearrangement. To incorporate this in the algorithm, as soon as 
the walker falls in a trap, we freeze that walk upto the next rearrange- 
ment. To ensure that a site once assigned as 'allowed', is not seen as a 
trap on a subsequent visit, we maintain a list of all allowed sites visited 
during r r , so that if these are revisited they are still assigned as 'allowed'. 

We average over 10 5 walks to get (r 2 ) for different time intervals 
t » r r .(r 2 ) vs. t are shown in fig 1, the curve is exactly like the results 
shown by ||. D(c, r r ) calculated from the average slope of (r 2 ) vs. t are 
plotted against r r for different c and shown in fig. 4. In fig. 5 we compare 
the simulated value of D with the value calculated from equation 12. 



4 Discussion 

Figure 5 shows that the calculated results are always lower than the 
value obtained from our simulations. The reason is probably that our 
simplified calculation ignores the distribution of distances travelled by 
different walkers and assumes an equal average distance covered by all 
the walkers. However the agreement is quite satisfactory inspite of this. 

The present model though closely related to the DBPM || and the 
finite trap model JT^| shows features distinct from either of them. We 



discuss here the points of similarity and dissimilarity between the models. 



Physically the difference between the finite trap model |T^] and the 
present one is only that here the traps are infinite. At this limit however, 
a qualitative difference enters which is the appearance of the percolation 
threshold. Now, in the quenched limit r r —>■ oo, when the trap positions 
are frozen, the diffusion coefficient is zero for nonzero c, so the threshold 
is at c = 0. The position of the threshold is at c = 0, whatever be 



the dimension of the system. So, the difference with the model |12 



is that there is a threshold. On the other hand the DBPM does have a 
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percolation threshold p c , which is at c = for 1-dimension, but for higher 
dimensions it is at the appropriate bond percolation threshold for that 
dimension. 

Our simulation results shown in figure 1 resemble the DBPM results 
for p < p c . So, the trapping effect is stronger here than in DBPM. This 
is expected since with sites as infinite traps, a random walker must fall 
in it given sufficient time, even if the trap concentration is very small. In 
the bond model 'forbidden' sites are inaccessible for the walker and may 
be avoided for c < p c . 

The finite trap model was shown to account successfully for the dy- 
namic disorder in the polymer electrolyte PEO — NH4CI4 @- There the 
crystalline regions in the polymer are assumed to be the sites with longer 
dwell time, and the amorphous regions the highly conducting sites with 
smaller dwell time. This model is suited for systems with stronger trap- 
ping, and also systems where reacting species are diffusing in a viscous 
medium. For example - formation for excimers in fluorescence exper- 
iments, where fluorescing monomers combine with each other to form 
an 'excited dimer' or excimer [HJ. The excimer has a certain lifetime 



after which it may dissociate, i.e. the walker is released from the trap. 
However, in this case dN oc c 2 rather than c in equation (5) if c is the 
monomer concentration. 

The primary interest of this model is that it is possible to calculate 
the dynamic disorder effect in a very simple way, compared to the DBPM 
which requires very involved mathematics, and the results are qualita- 
tively very similar. 
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Figure Captions 



Figure 1 : Plot of average squared distance versus time for c = .007 
from computer simulation. 

Figure 2 : Plot of fraction of walkers surviving versus time steps from 
simulation, equation (6) & equation (7) for c = .001 & c = .01 

Figure 3 : Plot of calculated diffusion coefficient versus trap concentra- 
tion for renewal time = 1, 10, 50, 250. 

Figure 4 : Plot of diffusion coefficient versus renewal time for c = .001, 
c = .005, c = .01 & c = .05 from simulation. 

Figure 5 : Comparison of diffusion coefficient versus trap concentration 
from simulation and theory [equation (12)] for renewal times = 250, 2000. 
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fig. 1 
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fig. 2 
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fig. 3 
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fig. 5 
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